Period fissioning and other instabilities of stressed elastic membranes 
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We study the shapes of elastic membranes under the simultaneous exertion of tensile and com- 
pressive forces when the translational symmetry along the tension direction is broken. We predict a 
multitude of novel morphological phases in various regimes of a 2-dimensional parameter space (e, v) 
that defines the relevant mechanical and geometrical conditions. Theses parameters are, respectively, 
the ratio between compression and tension, and the wavelength contrast along the tension direction. 
In particular, our theory associates the repetitive increase of pattern periodicity, recently observed 
on wrinkled membranes floating on liquid and subject to capillary forces, to the morphology in the 
regime (e<^.l,v^>l) where tension is dominant and the wavelength contrast is large. 



Thin membranes, such as paper sheets, tend to buckle 
when compressive forces are exerted on their boundaries. 
The origin of this familiar phenomenon, known as Eu- 
ler instability, is the large contrast between the energetic 
costs of straining and strain-free bending of thin elastic 
bodies. This feature is reflected in the different depen- 
dencies of the bending modulus B^Et 3 and stretching 
modulus Y^Et on the Young modulus E and thickness 
t of the membrane Wrinkling patterns, which often 
appear on supported membranes such as human skin or 
milk crusts, are characterized by a buckling scale Iq that 
can be much smaller than the membrane width W in 
the compression direction. Here, the distortion of the at- 
tached substrate gives rise to an energetic cost that is pro- 
portional to the amplitude Co of the bent shape, and the 
formation of a 1-dimcnsional (Id) periodic wrinkling pat- 
tern is induced by balancing restoring forces associated 
with this distortion and with bending resistivity of the 
membrane @. Near buckling (wrinkling) threshold, the 
ratio A — A/W between the displacement A and width 
W is small (see Fig. 1) , and the membrane attains a 
sinusoidal shape whose wavelength and amplitude are: 
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where o~ yy < is the compressive stress, and the inex- 
tensibility criterion (JT]b) reflects the conservation of the 
contour length of material lines along the compression 
direction y. This description corresponds to situations 
in which the exerted forces are purely compressive and 
enable the formation of a Id shape. This ideal picture 
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FIG. 1: Membrane geometry studied in this Letter. 



can be qualitatively modified when: (1) the translational 
symmetry along the orthogonal direction (x) is broken, 
or: (2) when tension a xx ~T is exerted on the mem- 
brane. The dramatic effect of symmetry breaking can be 
demonstrated by considering a bent membrane, whose 
amplitude (^(x,y) is forced to assume a sinusoidal shape 
of wavelength l e = 2o/(l + v) at the uncompressed edge 
(x = 0), where v is the wavelength contrast. As the in- 
cxtcnsibility criterion (JTJj) shows, the limit v—>oo cor- 
responds to a vanishing amplitude Ce - *• at the edge, 
a geometry similar to the "curtain problem" studied by 
Pomeau and Rica [|[ . As was shown in 0] , the transition 
from the vanishing amplitude at x~0 to the finite value 
Co, Eq. (jT]b) , is characterized by a cascade of sharp folds 
whose lengths may become infinitcsimally small as x^O. 
The important role of tension was recently noticed by 
Cerda and Mahadevan 0] , who showed that exerting ten- 
sion o~ xx =T on a membrane with length L and free edges 
at y=0,W implies a small compression \a yy \ ~ (t/L)T 
that may induce a periodic wrinkling pattern in the y di- 
rection. A useful number for characterizing tensile effects 
on wrinkling phenomena is the stress ratio e = \a yy \/o~ xx . 

In this Letter we classify the morphologies of clastic 
membranes for general values of the pair of parameters 
(e, z/) near buckling (wrinkling) threshold (A<1). Our 
theory reveals a surprisingly rich phase space, and in par- 
ticular predicts the existence of novel types of smooth cas- 
cades, qualitatively different from the hierarchy of sharp 
folds found in Q . Our interest in this problem was trig- 
gered by the repetitive increase of wrinkles periodicity 
in discrete steps, that was recently observed near the 
uncompressed edges (x = —L,0) of ultrathin membranes 
floating on liquid surface and compressed along y with 
A< 1 (||. As was realized in [3] , these membranes are 
under large tension o~ xx , and the amplitude £ e at their 
uncompressed edges is suppressed due to the high ener- 
getic cost of the liquid- vapor menisci (at x > 0, x < —L), 
induced by the bent membrane. The experimental con- 
ditions in [|| thus correspond to the asymptotic regime 
(e<Cl,^^> 1), which is the main focus of this Letter. 

Our geometry is depicted in Fig. 1. It is similar to the 
one studied in [3[ with three crucial differences: (i) The 
displacement A in the y direction is assumed constant 
(rather than constant compression a yy ~P). (ii) The 
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tension a xx = T in x is assumed constant (rather than 
constant length L), and (iii) The wavelength contrast v 
can be any positive number (rather than v— ►oo). In or- 
der to describe the wrinkling patterns of Q, we assume 
the membrane is floating on liquid, and include the en- 
ergetic cost of the lifted liquid mass of density p induced 
by the bent membrane. Our formalism and results are 
generalizable, however, to wrinkling problems in which 
gravity is replaced by another restoring force (e.g. due 
to an attached elastic substrate) or if bending is the only 
energetic cost (i.e. buckling). To quadratic order in A, 
the areal energy density of a shape ((x, y) is 0]: 

u= i(B(V 2 C) 2 +p. 9 C 2 +a(x)[(g) 2 -2A]+T(g) 2 ) (2) 

where cr(x), the compressive stress o~ yy that must be ex- 
erted at y=0, W in order to impose a constant A, appears 
as a Lagrange multiplier that forces inextensibility of con- 
tour lines parallel to the y direction. Let us briefly review 
the case of a free uncompressed edge y = 0, [j| . The 
energy is then minimized by the Id periodic pattern: 

(a) ax,y)=(sm(qy) ; (b) C = -V 7 ! (3) 

q 

with q = qo = (pg/B) 1 / 4 and a yy = <7o = —2-s/Bpg. The 
parameter e is hence defined as the ratio e = o~q/T. The 
energy of the shape ^ is the work ooAL, whereas the in- 
plane compression energy is hYA^VL. One thus obtains 
the threshold value A, n i„=2o'o/Y~, below which the mem- 
brane does not bend. Consider now the case in which 
the edge x = is forced to take an unstrained periodic 
shape ^p) with a wavenumber q e =(\-\-v)q$. We assume 
the length L is sufficiently large such that away from the 
forced edge at x=0 the membrane fully recovers its ener- 
getically favorable form (|3^) with q=qa- A natural guess 
for the shape is then the superposition: 

C{x,y) = (o(x) sin(q y) + sin(gij/) , (4) 
where: q 2 $ + q 2 ( 2 = 4A , (5) 
and: (a) lim d(x) = ; (b) lim ( (x) = , (6) 

x — >Xq x — >Xi 

where ([5]) is the inextensibility condition (recall that 
A <C f), and X o ,Xi = -L,0, respectively, and similarly 
qi=q e are introduced to simplify the forthcoming anal- 
ysis. Notice that the superposition ((4]) constitutes the 
most symmetric shape possible under the boundary con- 
ditions (BC), Eq. ©. Obviously, such a smooth shape is 
markedly different from the irregular one described in [3[ 
in the tensionless case (T=0) . This difference stems from 
the Gaussian curvature imposed in the membrane by a 
wavelength contrast v>Q 1 and the associated anharmonic 
energy density uq whose minimization (for T—0) gives 
rise to stress focusing at ridges and vertices which relieves 
the strain at all other areas of the distorted membrane 
. This principle underlies the emergence of sharp folds 
in the tensionless limit e— >0 0] . Assuming for the smooth 



shape §4$ d x (i^Q/l, with some typical length I, one ob- 
tains w G ~Fr 2 £iCM @,0|. ForT=S-Y^0 the situation 
is different since even regions free from Gaussian curva- 
ture are penalized by the term UT=T(d x () 2 ~Tr 2 J2i(i- 
A transition from an irregular shape (at e— >0) to a smooth 
one (jl]) is thus expected if ut > Ug- With Eq. © this 
inequality implies: A<5. Recalling the threshold condi- 
tion A > ctq/Y, we obtain a necessary condition for the 
existence of a smooth shape (U]): <5-e<A<<5^>e<f. 
We thus conclude that a transition between an irregular, 
sharply folded shape 0] and a smooth superposition ^ 
occurs at e SI ~0(1). A more careful analysis of Eulcr- 
Lagrange (EL) equations for the shape ^ reveals that 
e SI increases with the wavelength contrast v, and more- 
over e SI (y — > 0) > 0.5 0. Since for e, v— >oo the shape 
is characterized by a diverging number of generations of 
sharp folds 0], we conjecture: 

I. There exists a "branching" series {b n (e)}, such that 
for (e, v) with e>e SI (y) , b n (e)<i'<b n +i(e) , the morphological 
phase is I n , characterized by n generations of sharp folds. 

Let us focus now on the asymptotic regime (y, e <C 
e (is)), where a smooth shape described by Eq. (TJJ is 
expected. We transform now to the dimcnsionlcss set: 

- x _ y - C _ u _ a 

x=—\ y=—\ C= = — — ; u=— ; a=— ; <x,= veft'T, (7) 

where l T = s jTjpg and oo=l. Eq. |(2J) yields the EL Eqs.: 

M(a ,x)C - (I - 2ea 2 )^ + e%" = (8) 

M(ai, x)Ci - (I - 2ea 2 )Ci + <?Ci = (9) 
where: a 2 ( (x) 2 + a 2 Ci(x) 2 = 1 (10) 
and: M(a i: x) = a 4 + 1 + a(x)a 2 /e . (II) 

With the constraint ifTU]). Eqs. (JH® become a 4 th or- 
der nonlinear ODE for the function Co(x) (alternatively, 
£i(x)), which must be solved under the two BC ©. Ob- 
viously, two BC do not suffice to solve a 4 th order ODE. 
One may show, however Q , that in the regime 

e« 1 ; a a << 1/y/e , (12) 

the minimal energy profile that satisfy the BC © is de- 
termined, to leading order in e, by the 2 nd order ODE ob- 
tained from Eqs. (|8I9[) after neglecting the 4 th derivatives 
and the terms 2ea 2 Q . The BC © are thus sufficient for 
finding the minimal energy profile of the form ^ . Phys- 
ically, this means that force balance in the regime (Ti"2")) is 
dominated by tensile forces (— Q ) and by the restoring 
forces associated with variation of the pattern from its 
preferred wrinkling period (M(a,i, x)Q), whereas bending 
forces that result from variation along x (e 2 Q , 2ea 2 Q ) 
are negligible. Let us analyze now the asymptotic be- 
havior of Co (x) , £i (x) . In the limit x — > Xq (more generally, 
for x <$^. 1) we expect a(x)-^ — (a 2 1 +a ~ 2 )e=— 2e and from 
Eq. ([HD: M(a ,A > o)=0,M(a 1 ,A > o) = af + l-2a 2 >0. 
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FIG. 2: Computed patterns of wrinkled membranes (e<Cl), 
with wavelength contrast v=\ and a 1-strip shape (left), and 
with 1^=15 and a 4-strip series qo(l, 2, 4, 8, 16) (right). 

In this limit, linear analysis of Eq. @ implies: 

£ 1 = 9± sinh[JM( ail X )(x~X )} ; Co-^ 1 ^, (13) 
ai v a 

where Qi (oc e~ L ) is a constant that is determined by 
the solution of the nonlinear Eqs. ([8l fT0|) . The analogous 
analysis in the limit x — > X\ = relies on the assump- 
tion M(ao,Xi) < 0, that will be justified below. Linear 
analysis of Eq. ([8]) near x = X\ then yields: 

Co=— sm\J-M(a , Xi)(a-XJ]; Ci « 1 ~ , (14) 
a v oi 

where Zi is another constant. Some algebraic manipula- 
tions using Eqs. (|6I10I14I8I11[) yield: 

a 4 + 1 — a?{a\ + -k-) 
NUfli, a , Zi) = M (a ,Xi) = _ ° * ° ? (15) 

1 v ai / 

cr(Xi) = e[-(a 2 + ^) + ^i 2 M-(ai, a , Zi)](16) 

With Eq. (Tl~5jt . the asymptotes (|13I14|) . denoted hence as 
Cp(x), CP®! respectively, are fully described by the two 
unknown constants Qi,Z\. Although they are the lead- 
ing terms in asymptotic expansions whose radii of conver- 
gence are unknown (and may even vanish), we found an 
excellent agreement between the function obtained from 
Ci: CP through standard matching procedure and numer- 
ical solution of Eqs. (fSlfTTj) (in the regime (TT2"])) for all 
parameters ao,ai that we tested. Such matching analy- 
sis is based on equating Co(x) and (x), and their first 
and second derivatives at an unknown point x* < X\ . 
The three matching conditions translate into 3 algebraic 
equations for the unknowns x* , Z\,Q\, from which ex- 
pressions for Qi and Z\ as functions of ao,ai, valid in 
the limit e^Cl, are derived. Evaluation of the result- 
ing formulas Q proves that 0>i?i(ao=l, ai)>— 1 for all 
ai>l, and thus (see Eq. {HP ), confirming the assumption 
M(ao,-Xi) < 0. A characteristic shape of the form ((4]) 
with ai=2 is plotted in Fig. 2. The energetic cost U(v) 
(where v=a\— 1) of the shape (|U) with respect to the Id 
periodic wrinkling shape can be calculated (to leading 
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FIG. 3: Energetic costs U{v) (for e<Cl) of symmetric (1-strip), 
2-, and 3-strip shapes (obtained by minimizing over all series 
(l>ai>a2>l + ^). Dimensionless convention, Eq. iff}, is used. 

order in e) by substituting the obtained solutions in Eq. 
(|2|). We found the linear behavior: U~{y— 1), plotted as 
a gray solid line (logarithmic scale) in Fig. 3. 

The emergence of intermediate wavenumbers qo<qi<q e 
in the wrinkling patterns observed in Q inspired us to 
study the stability of the symmetric shape ([4]) with re- 
spect to "n-strip" patterns. A n-strip pattern is defined 
as n consecutive strips parallel to y. with a series of 
wavenumbers (qo < Qi - ■ ■ < (?n-i < q n — <Ze), and a series 
of borderlines (Xi, ■ ■ ■ ,X„_i) (with Xq=—L and X n =0), 
such that the profile in the j th strip is described by a su- 
perposition of the type (0]) with the change of subscripts: 

o -> j - 1, i -»• j ■ (17) 

As we now show, the shape in the j th strip, and in par- 
ticular its length Kj=Xj—Xj_i are determined iteratively 
from the shape in the (j— 1)*" strip by an energetic prin- 
ciple. Let us consider first the vicinity of the borderline 
X\ = -K\ in a 2-strip shape. We ask whether after the 
subscripts change (TP?)) .the shape can be described by the 
asymptotics (fT5|) (and Qz{x)=<S) at x— >K[, and by the 
asymptotics ([14")) (and Co(x) = 0) at x—*Kf, and a con- 
stant Q2- Such a shape implies a discontinuity of 
at x=Ki and hence the divergence of Q (x), for z=0, 2. 
These divergences, however, stem from neglecting the 4 th 
derivatives in Eq. I|8|). and are cured by the formation 
of a boundary layer of size e around x = K\ in which 
tensile and bending forces (— Q , e 2 Q , respectively) bal- 
ance each other. Similarly to the "take-off" line that is 
formed along a paper sheet when it is pushed into a nar- 
row ring, the energetic cost of this layer is negligible for 
e<l 0,11. Let us consider now Ci(x) near K\. Eq. ifTU]) 
implies Ci(-f i)=0 an d hence finite values of ( x (x) at both 
limits x — ► K^~, K^. If these values are different the 4" 1 
derivative e 2 ^ (x) diverges, leading to strong bending 
force that cannot be balanced by the finite tensile force 
(x), and whose existence thus leads to the emergence 
of a highly energetic region with localized Guassian cur- 
vature, similarly to ridges on crumpled papers @. The 
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"stitching condition": (x— > -KT)=Ci {%~ )■> amounts 
to the absence of such region. After manipulations simi- 
lar to those that led to Eq. (fT5|) this condition yields the 
asymptotics ([131) °f the second strip at x-^K± with: 



-(^fA£_(a 1 ,a ,Z 1 )Zl 



{18) 



M+(a 2 , ai,Q 2 ) = M(a 2 ,X 1 



i-{19) 



i + Qi(ff) 2 



The asymptotics of the second strip at cc— >A"2= are given 
by Eqs. (|14I17[) . The unknown constants are now K 2 = 
X 2 — X\ , Zq which can be found by following a matching 
procedure as described above, or by numerically solving 
the 2 nd order ODE obtained from Eqs. 1(51- 1101170 in the 
regime ((12|) on an unknown interval length K 2 and three 
BC: Ci(^2)=C2_(^i)=0 and C 2 (^ )=Q 2 \M^a 2 , ai , Q 2 ) . 
One thus obtains functional expressions for K 2 ,Z 2 in 
terms of ao,ai,a 2 . The construction of n-strip with a 
series (do = l<di<- • ■<a nr -i<a n =l+i>) proceeds iterativcly: 
The shape at the j strip is fully determined from the 
known shape in the (j — l) th strip by an identical pro- 
cedure to the one described above where the subscript 
change (jTTJ) is supplemented by 2 — ► j + 1. The energy 
of the n-strip shape is then computed from Eq. ([2]). 

In order to analyze the stability of the symmetric (1- 
strip) shape with respect to n-strip shapes, we calculated 
(for e<g;l) the minimal energies of 2— and 3— strip shapes 
for i>e(0, 20). Evaluation of similar plots for n>3 requires 
elaborate computations, since it involves (for every v) 
minimization in a high dimensional (d > 3) parameter 
space. We found that both 1— strip and 2— strip energies 
scale logarithmically with v , such that 1— strip is favorable 
for < v < /i (0), 2-strip is favorable for /i (0) < v < f 2 (0) , 
and 3— strip is favorable (over 1— and 2— strip) for v > 
/ 2 (0), where /i(0)«3.2, / 2 (0)«12. These observation 
leads to our second conjecture: 

II. There exists a "period-fissioning" series {/ n (e)} ; 
such that for (e,v) with e (y)<e<(L+i') 2 , /n(e)<^</r>fi(e), 
the phase is F n , characterized by n-strip shape. 

In Fig. 4 we depict the conjectured irregular ({I n }), 
symmetric (S) and period-fissioning ({F n }) phases in the 
2d parameter space (c.v). As will be shown elsewhere 
Q, the wrinkling patterns observed in [i| seem to be 
described by phases F 3 and F 4 . In order to understand 



the regime e>(l+^) 2 , consider Eqs. (1161170 . They predict 
that the compression |cr yj/ (a:) | in the j th strip increases 
roughly as q 2 , and hence imply that the local value of 
e (compression-to-tension ratio) increases and becomes 
0(1) for qj^qo/y/e- For wavelength contrasts v>t~ 1 / 2 
one may thus expect mixed phases (M nt k), characterized 
by n-strip that terminates at qj~qo / y/e, and followed by 
an irregular shape Ik that terminates at q e =(l+u)qo. 

The rich phase space depicted in Fig. 4 describes the 
morphology of stressed membranes under rather restric- 
tive conditions (A^;l, constant tension, and BC mapped 
onto a single number v). Future work will explore 
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FIG. 4: A conjectured phase diagram for stressed membranes 
(A<1,L>1). The axis correspond to the dimensionless 
parameters (e, v). The irregular (In), period-fissioning (Fn), 
and mixed (M n> k) morphologies are described in the text. 

whether morphologies beyond this regime are described 
by the conjectured phases or whether new phases emerge. 
Finally, the structure depicted in Fig. 4 raises the specu- 
lation that a nontrivial link exists between the irregular 
({I n }) and period-fissioning ({F n }) phases. Computa- 
tion of the conjectured series {b n (e), /n(e)} may reveal 
whether this impression is merely a superficial one or 
points to a deep duality between focusing and uniform 
distribution of stress in elastic membranes. 
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